Somatosensory evoked magnetic fields following stimulation of the ear in human 

T. Nihashi 1,2 , R. Kakigi 1 , O. Kawakami 1 , M. Hoshiyama 1,3 , K. Itomi 1 , 

H. Nakanishi 1 , Y. Kajita 2 , S. Inao 2 , and J.Yoshida 2 

department of Integrative Physiology, National Institute for Physiological Sciences, Okazaki, Japan 
2 Department of Neurosurgery, Nagoya University School of Medicine, Nagoya, Japan 
3 Department of Health of Sciences, Faculty of Medicine, Nagoya University, Nagoya, Japan 


1 Introduction 

Studies of somatosensory evoked magnetic fields 
(SEF) by averaging the signal of magneto¬ 
encephalography (MEG) have since revealed the 
homunculus structure of SI [1]: MEG is a non- 
invasive technique for the functional mapping of 
the human brain with high spatial and temporal 
resolution (in the order of mm and msec). This is 
accurate enough to estimate the location of the 
electrical source in the brain. 

The location of the face area in SI has been studied, 
mainly by stimulation of the lip, using electro- 
corticography in humans and monkeys. The 
location of the face area in SI has been also 
investigated by MEG [2,3], by positron emission 
tomography (PET) [4], and by functional magnetic 
resonance imaging (fMRI) [5] in humans. 
However, there have been no reports on stimulation 
of the ear. The objective of this study was to 
investigate the location of the ear area in SI. To our 
knowledge, this is the first systematic study of SEF 
following electrical stimulation of the ear [8]. 


intensity was approximately three times the sensory 
threshold (6-12 mA mean 10 mA). 
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Figure 1: a) Three points stimulated on the left ear. 
Black circle, black triangle, and black square 
indicates the stimulus point on the helix, lobulus, 
and tragus, respectively. 

b) A pair of silver ball electrodes. These stimulus 
electrodes were clipped onto the ear and fixed into 
position with a nose clip. 


2 Methods 

2.1 Subjects 

Thirteen healthy volunteers (8 males and 5 females; 
ranging from 25 to 46-year-old) were studied. 
Informed consent was obtained from all subjects 
prior to the study, which was first approved by the 
Ethics Committee at our Institute. 

2.2 Ear stimulation 

We stimulated 3 sites on the left ear: (1) the helix, 
(2) the lobulus and (3) the tragus (Figure la). A 
pair of silver ball electrodes were used for the 
stimulation (Figure lb). The stimulus electrodes 
were clipped onto the ear and fixed into position 
with a nose clip for swimming (Arena, Canada) or 
surgical tape. The electrical stimulus was a constant 
voltage square pulse wave delivered at a rate of 1 
Hz. We used a short stimulus duration (0.05 ms) to 
minimize electrical artifacts on SEF. The stimulus 


2.3 SEF recording 

The recordings were made in a magnetically 
shielded room. The subjects lay on the right side on 
a bed, and the SEF was measured with a 37-channel 
biomagnetometer (Magnes; Biomagnetic 
Technologies Inc., San Diego, CA). The center of 
the biomagnetometer was placed at the C4 position 
of the international 10-20 system, which covered 
both the neck and face areas of SI contralateral to 
the stimulation. We also placed the 
biomagnetometer in the ipsilateral hemisphere to 
the stimulated ear, but no good response could be 
recorded due to large artifacts. The MEG signals 
were filtered with a 0.1-100 Hz bandpass filter and 
digitized at a sampling rate of 2048 Hz. The 
analysis window was from 100 ms before to 300 
ms after the stimulation, and DC was offset using 
the pre-stimulus period as the baseline. Five 
hundred epochs were averaged in one trial, and at 
least two trials were recorded for each subject. The 




peak latency of each component was measured at 
the maximum value of the root mean square (RMS) 
of the magnetic field strength taken over the 37 
sensors. 

2.4 Single-dipole analysis 

A spherical model was fitted to the digitized head 
shape of each subject, and the location (x, y and z 
positions), orientation and amplitude of the best 
fitting single equivalent current dipole (ECD) were 
estimated at each time point. The origin of the 
head-based coordinated system was the midpoint 
between the preauricular points. The x-axis 
indicated the coronal plane with positive values in 
the anterior direction, and the y-axis the sagittal 
plane with positive values up. The correlation 
between the recorded measurements and the values 
expected from the ECD estimate was calculated, to 
evaluate how closely the measured values 
corresponded to the theoretical field generated by 
the model and the observed field. We estimated 
dipole fitting at a point with a correlation 
coefficient above 0.97. 

2.5 Double-dipole analysis 

We also used a brain electric source analysis 
(BESA) (NeuroScan, McLean, VA) computation of 
theoretical source generators in a 3-layer spherical 
head model. This method allows spatio-temporal 
modeling of multiple simultaneously activating 
sources over defined intervals. The location and 
orientation of the dipoles were calculated by an 
interactive least-square fit. 

When the data were read into BESA, the sampling 
rate was reduced to 850 Hz and a low-pass filter 
was applied automatically to 20 Hz. The analysis 
was performed as follows: In step 1, two dipoles 
were tentatively located in the neck and face areas 
of SI. The orientation of each dipole was calculated 
with the tentative locations. Then, in step 2, the 
locational constraint was released, and the locations 
of both dipoles were fitted, retaining the 
orientations obtained in step 1. Releasing the 
orientational constraint, the orientations of the 
dipoles were again calculated (step 3). The residual 
variance (%RV) indicated the percentage of data, 
which could not be explained by the model. The 
goodness-of-fit (GOF) was expressed in % as (100- 
%RV). We considered that when the GOF was 
larger than 90 %, the adaptation of the dipoles 
would be significant. We also tried changing the 
location and/or orientation of source 1 and/or 2 or 
to locate a new source, source 3, to improve the 
model. 


2.6 Statistical analysis 

We compared the peak latency of the early 
components of SEF among the stimulated sites. 
Statistical analysis was done by one-way factorial 
analysis of variance (ANOVA) followed by post- 
Hoe multi-comparison using Bonferroni-Dunn’s 
correction and P<0.05 was considered to be 
significant. 

2.7 MRI overlaying 

Magnetic resonance imaging (MRI) scans 
(Shimadzu Magnex 150 x T 1.5T, Kyoto, Japan) 
were obtained for all subjects. T1-weighted 
coronal, axial and sagittal images with continuous 
slices 1.5 mm in thickness were used for overlays 
with ECD sources detected by MEG. The same 
anatomical landmarks used to create the MEG 
head-based 3D coordinate system (the nasion and 
bilateral preauricular points) were visualized in the 
MRI images by affixing to these points high- 
contrast cod liver oil capsules (3 mm in diameter), 
the short relaxation time of which provides a high- 
intensity signal in T1-weighted images. 

3 Results 

3.1 SEF waveforms 

No clear and consistent component was identified 
in 6 out of 13 subjects, mainly due to large artifacts. 
Therefore, we analyzed results obtained from the 
other 7 subjects. One or 2 deflections were 
consistently identified within 60 ms in all subjects 
in the SEF waveforms following stimulation of 
each part of the ear (Figure 2). The peak latencies 
of the two components were approximately 20 ms 
(M20) and 40 ms (M40), respectively . However, 
there was a large inter-individual difference of 
waveforms as shown in Figure 3. We mainly 
analyzed the M20 component, since it was more 
clearly recognizable than the M40 in most subjects. 
However, in the subjects whose M20 was mixed 
with noise, the M40 component was used for the 
dipole analysis. The ECD direction of the M20 and 
M40 was posteriorly and anteriorly, respectively. 
Table 1 shows the peak latency of M20 and M40 
following stimulation of each part of the ear for all 
subjects. The latency difference was not significant 
for either M20 or M40. Only three out of 7 subjects 
showed a clear deflection at 80 ms (M80) after 
stimulation, no consistent response was identified 
later. 
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Figure 2: MEG waveforms in two representative 
subjects, subject 1 and subject 3, following 
stimulation of the helix, the lobulus and the tragus 
of the left ear. Three components, M20, M40 and 
M80, were identified. There was inter-individual 
difference in the waveforms. 


3.2 Analysis by single dipole model 

At first, we analyzed the ECD location for the M20 
or M40 component using a single dipole model. On 
stimulation of the helix, the ECDs of M20 or M40 
were estimated to be near the neck area of SI in all 
seven subjects. On stimulation of the lobulus, the 
ECDs were estimated in a similar area in four 
subjects. However, the ECD was estimated to be in 
the face area in one subject, and in the deep white 
matter in the other two subjects. On stimulation of 
the tragus, there was a large inter-individual 
difference among subjects in terms of the location 
of the ECDs. The ECDs were estimated to lie near 
the neck area in three subjects, near the hand area 
in two subjects, and in the deep white matter in two 
subjects. For three subjects who showed clear M80 
components, the ECDs were estimated to lie in the 
superior bank of the Sylvian fissure, corresponding 
to the secondary somatosensory cortex (SII). 

3.3 Analysis by multi-dipole model (BESA) 

Some subjects showed ECDs of M20/M40 in 
strange areas following the lobulus and tragus 
stimulation. In such cases, we used BESA to make 
a model with two dipoles locating in the neck and 
face areas of SI, source 1 and 2, respectively, since 
we speculated that the trigeminal and cervical 
nerves could be stimulated simultaneously. We 


obtained a reliably high GOF in all such condition. 
Figure 3 and 4 showed representative results for the 
BESA analysis, and the locations of the ECDs on 
2-D and 3-D MRI in one representative subject. 
The locations in the lobulus and the tragus 
stimulation in such conditions were estimated to be 
near the neck and the face area by the double¬ 
dipole model with a high GOF value. We also tried 
to make other 2-dipole or 3-dipole models placing 
sources at various sites, but the 2-dipole model 
described above was the most appropriate in all 
conditions. 
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Figure 3: Two-dipole model calculated by BESA for 
the responses to the electrical stimulation of the 
lobulus and tragus in subject 6. Source 1 and 2 
were located near the neck and face areas of the SI, 
respectively. Sources 1 and 2 corresponded to the 
M40 component and gray area showed the 
analyzing window for BESA. 



Figure 4: Two sources of M40 in subject 6 
overlapping on 3-D (b) MRI. 
















4 Discussion 

In neither human nor animals, had there been a 
systematic study of evoked responses following ear 
stimulation to represent the ear region in SI. To our 
knowledge, this is the first such study of SEF 
following electrical stimulation of the ear. 

One important finding was that two sources were 
activated in some conditions following lobulus and 
tragus stimulation. Using the single-dipole model, 
one source was estimated to be the middle of the 
two activities such as the hand area of SI or the 
deep white matter. Therefore, a two-dipole model 
using BESA was necessary to solve this problem. 
Dipoles were estimated to locate in the neck and 
face areas of SI. We speculated that the lobulus and 
tragus regions are innervated by the trigeminal and 
cervical nerves in some subjects. 

In developmental anatomy, the external ear 
develops from six mesenchymal proliferations at 
the dorsal ends of the first and second pharyngeal 
arches, surrounding the first pharyngeal cleft. These 
swellings (auricular hillocks) fuse and form the 
definitive auricle. The anterior three auricular 
hillocks mainly form the ventral part of the external 
ear and the posterior three auricular hillocks mainly 
form the dorsal part. The former is innervated by 
the auriculotemporaris nerve from V3 and the latter 
is innervated by the auricularis magnus nerve from 
the C2 and/or C3. Thus, it is possible that the 
lobulus and tragus regions receive sensory fibers 
from both the trigeminal and cervical nerves, 
although, in most adults, the cervical nerve 
innervation is dominant [6,7]. In the subjects with 
predominantly cervical innervation in lobulus and 
tragus areas, the ECD for M20/M40 was estimated 
to be located near the neck area. On the other hand, 
when the cervical nerve innervation was mixed 
with trigeminal nerve innervation, the “ear area” of 
SI showed a distribution both near the neck and 
face areas. In addition, there remains the possibility 
of a horizontal intracortical connection between the 
neck area and the face area. However, intracortical 
connections of somatosensory cortex were not 
confirmed in humans. 

In conclusion, the ear area of SI was located near 
the neck area, but the receptive fields of some parts 
of the ear, such as the lobulus and tragus, may be in 
both the neck and facial areas of SI. These findings 
suggested that the “ear area” of SI is variable unlike 
the other areas of SI, possibly because the ear is 
located on the border between the neck and face. 
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